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Abstract 



It has recently been argued that long range forces due to the exchange 
of massless neutrinos give rise to a very large self-energy in a dense, finite- 
ranged, weakly charged medium. Such an effect, if real, would destabilize a 
neutron star. To address this issue we have studied the related problem of 
a massless neutrino field in the presence of an external, static electroweak 
potential of finite range. To be precise, we have computed to one loop the 
exact vacuum energy for the case of a spherical square well potential of depth 
a and radius R. For small wells, the vacuum energy is reliably determined by 
a perturbative expansion in the external potential. For large wells, however, 
the perturbative expansion breaks down. A manifestation of this breakdown 
is that the vacuum carries a non-zero neutrino charge. The energy and neu- 
trino charge of the ground state are, to a good approximation for large wells, 
those of a neutrino condensate with chemical potential fi = a. Our results 
demonstrate explicitly that long-range forces due to the exchange of massless 
neutrinos do not threaten the stability of neutron stars. 
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I. INTRODUCTION 



If neutrinos are massless, the exchange of neutrino-antineutrino pairs gives rise to long 
range interactions between weakly charged particles, such as neutrons |1|J^ . On dimensional 
grounds, the potential between two neutrons separated by a distance r must be of the form 

.5 



(where is the Fermi constant), and thus falls off rapidly at large distances. As was first 
remarked by Feynman 0,^], an interesting effect may arise if there is a finite density of 
matter near the two neutrons. A 4-body process, for instance, can give rise to a potential 
that behaves like 

(1-1) 

if -R ^ r, where N is the total number of particles in the cloud of matter and R is the distance 
between the matter and the two test neutrons. Because of the large N"^ enhancement factor 
in ( |1 . 1| ) , the effective coupling GpN'^/R^ can be of order unity. In this case, far-away 
matter would give a large contribution to the two-body potential. Since this effect would be 
even more important for higher-order processes, one would then have to take into account 
diagrams of arbitrarily high order. 

More recently, Fischbach ||^ has studied the contribution to the energy of a neutron 
star due to the exchange of massless neutrinos. Large combinatoric enhancements, similar 
to those envisioned by Feynman, led him to conclude that a neutron star could not be 
gravitationally bound unless neutrinos are massive, rrii, > 0.4 eV.Q 

Abada et al have contended that Fischbach's calculation was incorrect. They have 
concluded that the neutrino contribution to the energy of a neutron star is essentially negligi- 
ble. While the arguments given by Abada et al were compelling, the calculations upon which 
they were based did not take into account the finite size of the star. As we will demonstrate 
in this paper, the conclusions reached in |^ were nevertheless essentially correct. 

We believe that we have completed the proof of the assertions made in 0. We have 
studied the ground state properties of a massless neutrino field propagating in a finite-range 
electroweak potential, like the one induced by a finite density of neutrons. To be precise, 
we have considered a spherical square well electroweak potential of depth a and radius R, 
and have computed to one loop the exact energy and neutrino number of the ground state. 
By exact we mean that our results are non-perturbative in the external potential. In order 
to carry out this program, we have used an approach due initially to Schwinger [|10|-Q For 



^ Smirnov and Vissani |p have argued that the long range forces can also be screened by the 
neutrino condensate that resides in a star 0. However, unlike screening at finite temperature, 
screening at finite density is not perfect. Fischbach, et al have argued, in response, that in a 
finite system a residual effect could still blow up a star (or make it collapse into a black hole). 

^Similar calculations are found in non-perturbative computations of quantum corrections to soli- 
tons [pl-|l3|| and in studies of QED in strong fields 
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technical reasons, the potentials which we consider are of relatively modest strength and 
range compared to those of a realistic neutron star. As we shall argue, our calculation is 
nevertheless sufficient to demonstrate that the neutrino ground state energy in a neutron 
star does not vary wildly as the radius or density of the star is increased, but is in fact very 
well-behaved. 

The paper is organized as follows. In Section H we introduce the reader to Fischbach's 
results and review the critique advanced by Abada, et al. We then make the case for our 
approach to the question. We show that we have a well-defined problem to which we can 
give a definitive solution. 

Section |T| is essentially independent of the specific application to the neutron star prob- 
lem and, we believe, is of interest in its own right. We first solve the Weyl equation in the 
presence of a spherical square well potential. Because the fermions are massless, there are 
no bound states but only scattering solutions. From the expression for the phase shifts, 
we compute the properly renormalized energy of the neutrino ground state in the presence 
of the potential. For small wells, aR <^ vr, we show that the energy of the ground state 
is perturbative in the potential. As the size of the potential increases, the perturbative 
expansion breaks down. A symptom of this breakdown is that for aR > ir the ground state 
carries a non-zero neutrino charge. We show that, to a very good approximation for large 
wells, the ground state energy and neutrino number are those of a neutrino condensate with 
chemical potential fi = a, a result anticipated by Abada et al. That this result also holds 
for a neutron star follows trivially. 

We emphasize that we have studied precisely the same problem as Fischbach did. Our 
approach, however, is quite different. We conclude that long range forces due to the exchange 
of massless neutrinos do not destablilize a star. There is indeed a shift in vacuum energy, 
but - maybe unfortunately - it is too tiny to be of any consequence regarding the fate of a 
neutron star. 

Related lines of argument can be found in the recent literature. In ||T3[, a 1+1 dimensional 
star is studied as a toy model. In this case, the vacuum energy can be calculated exactly and 
happens to vanish. It is also shown that the presence of a boundary implies the existence 
of a neutrino condensate. This point is also emphasized in [0. Both papers stress that a 
complete 3 + 1 dimensional calculation for a system with a boundary is necessary in order 
to settle the issue definitively. 



II. THE NEUTRINO GROUND STATE IN A NEUTRON STAR 

A. Preliminaries 

Before embarking on our main calculation, we introduce some useful concepts and for- 
mulas. The stated goal in Ref. is to calculate the shift in vacuum energy due to neutrino 
exchange in a neutron star: 

Ty=(0|i/|6)-(0|i7o|0). (2.1) 

Here H denotes the Hamiltonian in the medium and \0) refers to the true ground state of the 
system. The second term is the usual matter-free vacuum energy. For fermions, Eq. ( p.l|) 
amounts to 
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^ U>o j<o ) 



(2.2) 



where the indices i refer to the positive and negative energy levels of the neutrinos. A first 
step toward the actual computation of ( p.l|) is the following expression due to Schwinger [jlO 



W = - J d^x — ti [Spix, x') - SUx, x')] (2.3) 

where Sf{x,x') is the 2x2 Feynman propagator for the massless left-handed neutrino field 
in the presence of the finite neutron density and Sp{x,x') is that of the free field. In ( p. 3D , 
the limit x' — x is taken symmetrically and averaged. 

Let us now drop any explicit reference to the neutrons of the star and consider instead 
the following effective Lagrangian for the neutrino field 

£eff=?^L[^^+«7°] (2.4) 

where ipi = Lip, L = {1 — 75)/2, and where 

a{x) = GFpn{x)/V2 ~ 20 eV 

is the electroweak potential induced by the finite neutron density (p„ ~ 0.4 fm~^ in a 
typical neutron star). The static potential a > is attractive for neutrinos and repulsive 
for antineutrinos. The effective Lagrangian of (|2.4|) has been derived in various ways from 
the underlying electroweak theory [|1^-[19| and applied to studies of the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect It provides us with a most convenient tool to study the 

coherent behavior of low-energy neutrinos in a neutron star, on scales large compared to the 
size of a neutron r„ ~ 1 fm (200 MeV)~^.0 

Since the potential is static, it is convenient to introduce 

SP {x, x') = — f duj e-^^(*-*') G(o) (x, x'; u) (2.5) 



27r Jc - -'^'^ 

with C representing the usual Feynman contour. The resolvents G{x, x'; uj) and Gq{x, x'; uj) 
satisfy 

{iu - i/(o)) G(o) (f, x'; u) = -6%x - x') (2.6) 



with the Hamiltonian from ( |2.4| ) 

H = -a + ia-V,c] Ho = ia-V^. (2.7) 
The Schwinger formula, (|2.3|), becomes 



^ The fact that the effective theory, as in (2.4), is anomalous should not worry us. The underlying 
theory (the Standard Model) is perfectly well defined. Besides, for a static potential in 3 + 1 
dimensions, there are no anomalous efi'ects because E ■ B = 0. 
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W=^ I duuTi^[G{uj)-G,{uj)] 

ZTT JC 



27r Jc 



duj uj Tr^ 



CO — H CO — Ho 



where the trace is over spinorial and configuration space indices.^ Integrating 
and expanding the logarithm in powers of the external potential finally gives 



W 



/ dto Trx log 


co-H- 


.uj - Ho. 


Jc 



1 °° 1 /■ 



(2.8) 

by parts 

(2.9) 
(2.10) 



Even with the free particle vacuum energy subtracted, Eq. ( |2.10| ) is still a formal, 
ultraviolet-divergent quantity. The culprit is the first term,0 k = 2, which is the famil- 
iar vacuum polarization diagram. This diagram is superficially quadratically divergent and 
has to be renormalized. However, because of gauge invariance, the actual ultraviolet di- 



vergence is milder. We will discuss this point in more detail in Section [lTl| . Let us simply 
mention here that 



strictly. As R increases, the k = 2 term gives a vanishing contribution to the energy per 
unit volume of the star.Q 

The remaining contributions, k > 4, are ultraviolet convergent because the interaction 
term in ( ^.4[ ) is renormaiizable. On the other hand, because the neutrinos are massless, 
these terms have a strong infrared dependence 



k R 



(2.11) 



For a neutron star, a ~ 20 eV and i? ~ 10 km, so that aR ~ 10^^. With such an expansion 
parameter, the perturbative expansion of (2J.) is doomed to diverge and Eq. ( p. 10 ) must be 
resummed. 



The rule for closing the Feynman contour C in the complex uj plane in Eq. (2.8) is to take the 



average of the integral around the positive and negative energy cuts, so as to recover (2.2). 



^Only the k — even terms contribute for a static external potential. This is most easily seen in 
the Dirac representation of ( |2.4| ) because then the coupling to the potential involves the 75 matrix 
and the trace over an odd number of vertices is proportional to the Levi-Civita tensor. For a static 
potential, there are not enough indices to be contracted with the Levi-Civita tensor and the odd 
terms vanish identically. 

^This statement is correct for a smooth effective potential, a. In a realistic star, the coarse grained 
structure at the neutron scale leads to a vacuum polarization contribution which scales like R^. 
See also the "Note added in proof." 
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For comparison, the estimate of the vacuum energy (|2.1|) given in P] is a large but finite 



sum, 

N 



W = J2W^''\ (2.12) 

k=2 

where ~ 10^'^ is the number of neutrons in the star. Each term is taken to be given by the 
average k-hodj potential energy multiplied by the number of ways of choosing k neutrons 
from the grand total of A^. For k large, but much less than A^, 



Since aR ~ 10^^, it is not surprising that after only a few terms the sum in Eq. (2.12) leads 



to an estimated energy which exceeds the mass of the neutron star. For future reference, let 
us quote the following estimate obtained in Ref . . For fixed neutron density, corresponding 
to a ~ 20 eV, a subsystem of radius i? ~ 2 ■ 10^^ cm contains as much energy as the total 
mass of a neutron star with a radius of 10 km. Note also that aR = 20 eV x 2 ■ 10~^ cm ^ 20. 

Although similar, the estimates ( p.ll|) and ( p.l3| ) differ crucially in the details. As 
emphasized by Abada, et al 0, for fixed fc, there are 0[{k — 1)!] terms missing in ( p. 131) . 



These are associated with multiple rescattering over the same neutrons. These processes 
also turn the finite sum ( |2.12|) into an infinite series, as in ( |2.10|) . For large aR, the series is 
diverging and must be resummed. 



B. The scattering problem 

Consider the expression given in (|2.2|) for the vacuum energy. We can count the energy 



levels of the Hamiltonian (|2.7] ) by putting the system in a large box Let us specify each 
eigenstate by its energy eigenvalue E and by a set of quantum numbers k. For simplicity, 
consider a spherically symmetric potential in a large spherical box of radius -Rbox- (The 
parameter k is then related to the total angular momentum.) Because the neutrinos are 
massless, the spectrum of ( p.7|) has no bound states but only scattering solutions, 

V'L« ~ sin(/cr + 0«), (2.14) 

where k is the radial momentum. At the boundary of the box of radius -Rbox impose, for 
instance, that 



fc(o)i?box + = rivr n = 0,1,2,... (2.15) 



where the index refers to the free modes. From ( |2.15| ), the shift of energy level for fixed k 
is then given by 

„ „ dE , , , . 1 dE , n N dE Ok 

E-E.^ _(, -k,,)^-—-^ i,. - fit f^-"' 

where 5^ is the phase shift between the free and scattering eigenstates. Given that the 
number of modes per unit energy is 
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we finally get 



dn_l dk 

dE - ^^-^^^ 



dn{E - E^) = --5^dE. (2.18) 

TT 



Hence, Eq. ( p.2|) becomes |jTO| 



^ = 7^E/ dE[5^{E) + 5^{~E)]. (2.19) 
27r "T^ JO 



Of course, just as in ( |2.10| ), the expression Eq. ( p.l9| ) is a formal, ultraviolet-divergent 
quantity that must be renormalized. 



Our strategy to address this issue is the following. (Details may be found in Section |T|, 
where we carry out an actual calculation.) Consider the expansion of the vacuum energy 
in the external potential of Eq. (|2.10|) . As noted above, only the term of second order 



in the external potential is ultraviolet divergent. The evaluation of this term involves the 
vacuum polarization tensor II^i,. Its divergent part can be absorbed using the standard 
renormalization procedure. The terms with > 4 in the expansion ( p. 101 ) are ultraviolet- 



convergent. If aR 3> 1 (with R the effective range of the potential), however, the terms 
with /c > 4 in ( [2.101 ) must be resummed in order to get a sensible result. Using Eq. ( ^.lOp , 
this seems like a hopeless task. 

Consider then the Taylor expansion of the phase shifts 5^ in powers of the external 
potential aR. Only the first and second Born terms will lead to ultraviolet divergences in 
the expression of the vacuum energy ( |2.19 ). The first Born term is odd in E, and drops from 



the calculation if we integrate symmetrically over positive and negative energy phase shifts. 
If we subtract the second Born approximation from (p.l9|) , the resulting, "renormalized 
vacuum energy" is free of ultraviolet divergences. The resulting quantity is equal to the 
resummed series of Eq. (|2.10|) with k > 4 and is precisely the quantity which has been 
estimated in Ref. [Q. 

In Section |III| we compute the renormalized vacuum energy for a spherical square well 
potential with depth a and radius R using Eq. ( [2.191 ). We numerically integrate the phase 



shifts over energy and subtract the second Born term. In principle we are required to evaluate 
all of the terms in (|2.19|) . In practice, it is sufficient to investigate the convergence of the 



series in order to obtain a numerical estimate for the vacuum energy. As we will see, for 
large aR the number of terms to be calculated grows linearly with aR. Of course, we cannot 
expect to carry out such a calculation for an actual neutron star, for which aR ~ 10^^. This 
is not necessary, however. Because the neutrinos are massless, there is just one dimensionless 
parameter, which is the effective expansion parameter aR. The question is then how the 
ground state energy changes from the domain where the perturbative expansion converges, 
aR -C 0{1), to the domain where the expansion presumably breaks down because of the 
infrared divergence, aR ^ 0{1). 

We have computed the renormalized vacuum energy for a number of points in the range 
< aR < 100, including the point aR = 20. According to Fischbach, et al 0, the 
energy corresponding to this point should be on the order of the mass of a neutron star. 
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W ~ 10^° kg ~ 10^^ eV. By way of comparison, we find W ~ —2.3 keV. The energy is well- 
behaved over the whole range of our calculation. Remarkably, the energy per unit volume 
of the potential plotted as a function of aR exhibits a crossover at aR = vr. For aR > vr, 
the ground state of the system carries a non-zero neutrino charge. 
The neutrino charge of the vacuum in the potential a is defined as 




; J d^xti [Sf{x,x') - S^p{x,x')]^,^^ (2.20) 

(2.21) 



using the notation introduced in (|2.2|) . Using Eqs. (|2.15|) - (|2.18|) gives (|2.2CI| ) in terms of the 
phase shifts, 

,._i^r,sM_^|. (2.22) 

^ 2tt^Jo \ dE dE j ^ ' 

Just as was the case with the formal expressions for the ground state energy, the above ex- 
pression for the vacuum charge will need to be renormalized. In this case this is accomplished 
by subtracting out the first Born term. 

The existence of a neutrino condensate in the ground state is to be expected . Because 
the neutrinos are massless, neutrino-antineutrino pairs are a priori easily produced by any 
non-zero gradient of the potential. Perhaps surprisingly, we have found that in our model it 
takes a finite value, aR = tt, for this to happen. For aR > vr, the ground state of our model 
contains a net neutrino charge. Note that the coupling term in ( p.4|) is precisely analogous 
to a constraint on the neutrino number with chemical potential n = a |]T8[. In the large 



volume limit, we expect the charge of the vacuum to be well approximated by that of a 
condensate with volume V = inR^/S [^, 

Similarly, using the dispersion relation uj{k) = k — a [which follows from ( p.4|) in the case of 
a homogeneous potential] , the renormalized ground state energy should approach [] 

^cond = VJ ^(fc - a)9ik -a) = (2.24) 

That this is indeed the limiting behaviour of the renormalized vacuum energy for large 
potentials is demonstrated in Section pTl] . 



^In a neutron star, the induced neutrino charge density is extremely small, ~ x2- 10~^^fm~^ 
(for a ~ 20 eV, and with the number of massless neutrino species), compared to the neutron 
density, ~ 0.4 fm~^. Consequently the back-reaction of the neutrino condensate on the neutron- 
induced effective potential is totally negligible. Also, because the low energy neutrinos of the 
condensate interact very weakly with each other, higher orders effects are presumably very small. 
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III. THE NEUTRINO GROUND STATE IN A SQUARE WELL POTENTIAL 



A. Derivation of the phase shift 

Both the charge and the energy of the neutrino ground state in the presence of a back- 
ground potential may be expressed in terms of the scattering phase shifts. Our starting point, 
then, is the equation of motion following from the effective Lagrangian given in Eq. (|2.4|), 



uj-ia-V + a{r)\'ipL = 0, (3.1) 

where ipL is a two-component spinor, and where the potential is taken to be spherically 
symmetric. The energy, uj, may be positive or negative, with positive values corresponding 
to particle solutions and negative values corresponding to anti-particle solutions. 



In order to solve for the scattering solutions we proceed, in the usual way by 
decomposing them into angular momentum eigenstates. 



where 




^L,. = fjmir)n\i! + g,^ir)n^-J, (3.2) 



^ = J - 2 (3-3) 

l'=J + l>0, (3.4) 

and where we use the standard notation for the spherical harmonics. The two-component 
spinors fl^j^ are eigenstates of and Jz, with eigenvalues j{j + 1) (j=|, |, . . .) and m, 
respectively. They are also eigenstates of L^, with orbital angular momentum l=j — | for 

flj^ and l'=j + I for fi^"''. In the following we will label the solutions by /= 0, 1,2, . . . 
instead of by j, although it is to be understood that the solutions themselves are eigenstates 
of the total, not the orbital, angular momentum. With the above normalization for the 
angular momentum eigenstates we have the useful relation 

a ■ f = ^^jm- (3-5) 

Inserting ( |3.2| ) into (|3.1|), we then obtain a set of coupled first-order differential equations in 
// and gi, 

^ --fi = + a{r)]gi (3.6) 
ar r 

^ + ^-^gi = -i[uJ + a{r)]fi. (3.7) 
dr r 

As noted above, we will use a spherical square well potential in order to simplify our 
calculations. This choice has the benefit that the expression for the phase shift may be 
written out explicitly in terms of simple functions. We thus take 
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a{r) = a9{R-r), r > 0. (3.8) 



The solutions of Eqs. (|3.6|) and (|3.71 ) are then simply given by spherical Bessel functions. 
Requiring that the solutions be regular at the origin, and dropping an over-all normalization 
factor, we find 

ji{\LU + a\r), 0<r<R 

fi = { (3.9) 

B [ji{\Lu\r) cosSi + ni{\uj\r) sin 5;] , r > R 

—ie{ijj + a)ii^i{\uj + a|r), < r < 

91 = { (3.10) 
-ie{uj)B [ji+i{\uj\r) cosSi + ni+i{\uj\r) sin 5;] , r > R, 

where B and 6i are fixed by matching /; and gi at r=R (note that the derivatives of // and 
Qi are in general discontinuous across the boundary) and where 

e{x) = ^. (3.11) 
\x\ 

Our conventions for the spherical Bessel functions are as follows: 

Jiip) = (^) ' Ji^ip), niip) = i-iy (^^) ' J_,_. (p). (3.12) 

Performing the matching at the boundary, we finally arrive at the desired expression for the 
phase shift 

, X ^ g(^)j;(k + (^\R)ji+ii\^\R) - + a)ji{\^\R)ji+i{\^ + a\R) /o 1o^ 

'''' '^""^ £(a; + a)nK|a;|i?)jz+i(|a; + a|i?)-£(a;)jK|^ + a|i?)nm(kl^)' ^' ' 

We may now use this expression in order to calculate the energy and charge of the vacuum. 



B. Energy and charge of the vacuum 



The formal expressions relating the energy and charge of the vacuum to the scattering 
phase shifts are given in ( |2.19D and ( p. 221 ). The spherical symmetry in our problem implies 
that the sum over n becomes a sum over /. We may thus rewrite these expressions as 



1 r'^ 
W=—J2{2l + 2) du [6iiuj) + 6ii-u)] 
2tt Jo 



d 



.^-^E(2/ + 2)f 



(3.14) 
(3.15) 



The factor {21 + 2) = (2j + 1) is the degeneracy factor for a given energy u) and total angular 
momentum j. 

As we have discussed in Section ||, the above formal expressions are in need of renor- 
malization. Since our model is renormalizable, the divergences are confined to the first few 
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Bom terms. The procedure to be followed in handling these divergences is then as follows: 
(i) Taylor-expand the formal expressions for the energy and the charge in order to isolate 
the divergences; (ii) subtract out the divergent terms and "resum" the Taylor expansion; 
(iii) renormalize the divergent terms in the usual way and add them back in. The resulting 
renormalized expressions are then finite and correspond to the actual energy and charge of 
the ground state of the system. 

The first step in renormalizing the energy and the charge is to perform a Taylor expansion 
of the phase shift in aR. Let us then define 



5i{u- aR) = {aR)6l^\uj) + {aRf5\^\uj) + {aRf5f'\uj) + 



(3.16) 



where 



1 Ff^ 



(3.17) 



The explicit expressions for the first two Born terms are given by 



(3.18) 



nimm + nz+i(/?)j,+i(/?) - - {Ij^m+m + {l + 2)ji^,{(3)nim 



X 



2(/ + i: 



(3.19) 



where we have defined (3=\uj\R. It is easy to convince oneself that the even (odd) Born 
terms are even (odd) in uj. 

The formal expression for the energy is rendered finite by subtracting out the second 
Born term. (Since the first Born term is odd in cj, it drops out of the calculation if we 
integrate symmetrically over positive and negative energies.p|) To compute the second Born 



term, we may use the expression for the k = 2 term in the expansion (|2.10|) . Using the 
standard prescription, we may then regularize and renormalize this term. The renormalized 
expression for the energy is then given by 



(3.20) 



where 



iy(4+) = _ y(^2l + 2) duj \5i{uj) + 5i{-u) - 2{aRf5\^\u) 
27r r-^ JO L 



(3.21) 



^Note that all terms in the Born expansion which are odd in a are also odd in oj and thus do not 
contribute to the energy. See the footnote following Eq. ( |2.1C1| ). 
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and where W^^^ denotes the renormahzed vacuum polarization contribution. corre- 
sponds to the resummed terms in the expansion ( p. 101 ) with k > A. 

The formal expression for the charge in (|3.15|) is divergent because the phase shifts tend 
to a constant as a; — > oo 



roo 

lim Si{uj) = ± / a{r)dr = ±aR. (3.22) 

L0^±QO Jq 

This divergence is removed by subtracting out the first Born term, which, as may easily be 



verified from (|3.18|) , satisfies 



lim 6l^\uj) = ±1, 6l^\0) = 0. (3.23) 

Once regularized, this term gives a vanishing contribution to the renormahzed charge (the 
charge is actually zero to all orders in perturbation theory), so that the final expression for 
the renormahzed charge is given by 

1 oo 

^ren = — E(2^ + 2) fe(0+) " HOl] ■ (3-24) 

2^ 1=0 ^ ^ 

The problem of calculating the charge of the vacuum then simply reduces to that of evalu- 
ating the particle and anti-particle phase shifts at the origin. 

It is useful to note that for small aR the Taylor expansion of the phase shift given in 
( p.l6| ) may be used to obtain a perturbative expansion for W'^'^~^\ For future reference, let 
us write this perturbative expansion as follows: 

-1 oo 

Vr(^+) = — - J2 («^)"E(2^ + 2)^z"\ (small ai?) (3.25) 

2^-^ n=4,6,8,... 1=0 

where Ii^^ denotes the integral over oo of the n^^ Born term, 

li""^ = 2 r d{ujR)Sl''\uj). (3.26) 



This expansion is expected to be reliable for small aR. For large aR, of course, the pertur- 
bative expansion breaks down. 

The main goal in the remainder of this paper will be to study the behaviour of the 
renormahzed energy and charge as a function of aR. As we shall see, for fixed a both 
of these quantities scale like the volume of the potential region as the volume gets large. 
Before studying these quantities in detail, however, let us first consider the contribution 
to the total energy due to the vacuum polarization, W^^^. Since W^^l scales roughly like 
the surface area of the potential region (for fixed a), its contribution to the total energy 
is eventually overwhelmed by that due to W^^^\ Nevertheless, the handling of this term 
is somewhat subtle, so it is important to consider it carefully. Once we have studied W^^^, 
we will dispense with it altogether and consider only the contribution due to when 
discussing the renormahzed energy. 
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C. Vacuum polarization 



We now turn to the calculation of W^4n- Consider the k = 2 term in the expansion ( ^.lUp , 



Using 



leads to 



ly(^) = f du f d^xid^X2 tr [a{xi)Go{xi,X2)a{x2)Go{x2,Xi)] . (3.27) 
Am J J 



^3 h, pih-x 



"^'^ a{xi)a{x2)Tioo{p)e'^-^^'-^^\ (3.29) 




<i ^ - 7 + log47r + ^ - log (^-|^^ I (3.30) 

is the familiar one-loop vacuum polarization tensor, which we have calculated using dimen- 
sional regularization. In Eq. ( p.30| ), p is the 4-momentum in Minkowski space and ^ is the 
renormalization scale. The divergent term, oc 1/e, (plus an arbitrary constant term) can 
be absorbed as usual by redefining the parameters of the theory, a — > a(/i). Note that, 
because the polarization tensor is proportional to p^, the renormalized vacuum polarization 
term vanishes if the potential is homogeneous 0. Introducing 



a{\p\)=aj d\e'^^-'e{R-r) = ^^j,{\j^R)^ (3.31) 
with ji{x) = sin(x)/x^ — cos(x)/x, gives finally 



As it stands, the renormalized expression in ( p.32| ) is still a linearly divergent quantity. We 
have chosen to work with a spherical square well potential and, because of the presence of 
the sharp edge, arbitrarily high radial momenta contribute to the integral. Any smoother 
potential would lead to a finite result. At the level of the effective theory, there is a natural 
ultimate cut-off provided by the size of the neutron, A ~ 0.2 GeV. (This also suggests the 
choice fi = A.) Inserting the cut-off, we get that the vacuum polarization term ( |3.32| ) scales 



roughly like the area of the potential, VT^n ~ A. Consequently, W^^^/R^ contributes a 
vanishing amount to the energy density as the size of the system increases. Let us emphasize 
that the remainder of our calculation is completely insensitive to the presence of this cut-off 
since the higher order terms in ( |2.10| ) are ultra-violet convergent. 
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D. Numerical evaluation of the renormalized energy and charge 

We turn now to a numerical evaluation of the renormalized energy and charge. From 
now on we will ignore the vacuum polarization contribution to the energy. We will thus 
loosely refer to as the renormalized energy. 



1. Small values of aR 

Let us begin our study of the charge and energy of the ground state by considering rather 
modest values of aR. Fig. ^(a) shows plots of the particle and antiparticle phase shifts (solid 
curves) as functions of u for / = in the case that aR = 1. Also shown are the sums of the 
first and second Born approximations to the phase shifts (dashed curves). Note that the 
particle (antiparticle) phase shifts approach +{—)aR as — ^ oo. Furthermore, the phase 
shifts go to zero at the origin so that, according to Eq. ( p.24|) , the vacuum charge is zero.0 
The plots for higher values of / are qualitatively similar. 

Fig. |l](b) shows a plot - again for aR = 1 and / = - of the sum of the particle and anti- 
particle phase shifts with the second Born term subtracted. This is the quantity which must 
be integrated (and summed over I) in order to obtain the renormalized energy in Eq. ( p. 211) . 



For small enough values of aR, one expects the energy to be well-described by the leading 
term in the perturbative expansion defined in (p.25| ). Our numerical results show that this 
is indeed the case for aR < 1/2, which actually provides a non-trivial check on our results. 
The advantage of using the perturbative expansion to evaluate the renormalized energy for 
small aR is that the integrals over u may be done exactly in that case,|3 whereas those 
for the exact expression, ( p.21| ), must be done numerically. Furthermore, in the case of the 
exact expression, the numerical integrations over u become increasingly difficult as aR — >■ 0. 

The leading term in the perturbative expansion is at fourth order in aR. Fig. |^ shows 
a log- log plot of (—1) X (2/ -|- 2) X j/^'' as a function of I. (We have not plotted the point 



^ It may at first seem surprising that the charge vanishes for aR = 1. After all, the neutrinos are 
massless and it would seem that any gradient of the potential could produce neutrino-antineutrino 
pairs. However, this argument neglects the effect of the neutrino's spin. To see this, consider 



the system of first-order differential equations for fi and gi, (3^) and (p7?|). Writing these as two 
decoupled second-order Schrodinger equations shows that the gi component always has a repulsive 
centrifugal barrier, 

-FT + ~^ 2 — + ^ + a) /z = 

ar^ r dr 

d^gi 2dgi (/ + !)(/ + 2) ^ , ^ .2 n 

-TT + 2 91 + (w + a) 5z = 

ar^ r dr 

Thus, even for the / = mode, it takes a finite aR to create a charge in the vacuum. 

^'^We do not include the explicit form of 5\^\uj) here because it is rather unwieldy. As was the 
case for 5^\ijj), however, the expression for 5\'^\uo) may be expressed in terms of spherical Bessel 
functions. Thus the integrals may be done analytically (albeit with the aid of a computer). 
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corresponding to / = for obvious reasons.) The sum over / of this quantity is directly 
proportional to the renormalized energy in the perturbative region. As is clear from the 
figure, this quantity has roughly a power-law fall-off as / — >■ oo. Using the slope of the 
curve in Fig. ^, we find that {21 + 2) x //^-' oc for large /. This allows us to obtain an 
approximate answer for the sum over /, 

oo 

^(2/ + 2)4^^ = -0.00721 ± 0.00005, (3.33) 
/=o 

where the quoted uncertainty is a rather conservative estimate of the error incurred in the 
extrapolation to large /. We are thus able to compute the leading perturbative contribution 
to VT*^^^-*, which is given by 

= _o.00115(ai?)Vi?+... (3.34) 



2. Larger values of aR 

As aR is increased, one finds that there is a critical point beyond which it becomes 
energetically favourable for the neutrinos to condense; thus, for aR>aR\^^-^^, the ground 
state of the system contains a neutrino condensate. The presence of the condensate has a 
very interesting effect on the renormalized energy. Before calculating the energy, however, 
let us first examine the phase shifts for large aR. Fig. ^ shows plots of the particle and 
antiparticle phase shift for aR = 10, with / = 0,2,4 and 6. These plots are qualitatively 
very different from those for small aR [c.f. Fig. 0(a)] in that the phase shift in the strong- 
potential case exhibits resonances. Furthermore, for / less than some critical value the phase 
shifts do not vanish at the origin, as they do in the weak-potential case.0 Thus, according 
to Eq. ( p.24| ), the ground state is charged: 

QreniaR = 10) = ^ = 56. (3.35) 

1=0 

Note that the sum over / for the charge always truncates at some lma.x- For aR = 10, 
/ — 5 

The charge of the ground state as a function of aR is intimately tied to the resonance 
structure of the particle and antiparticle scattering solutions. It is well-known in non- 
relativistic scattering theory that as a potential well is deepened, the resonances of the scat- 
tering solutions migrate to the origin of the complex fc-plane and, at some critical strength 
of the potential, join the imaginary axis and become true bound states. As the potential 
is made deeper and deeper, the bound states move to lower and lower energies (i.e., they 
become more and more tightly bound). The situation is slightly different in the relativistic 
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Note, however, that Levinson's theorem is still satisfied: since there are no bound states, the 



sum of the particle and anti-particle phase shifts at the origin is zero [23|. 
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scattering theory for fermions. Consider first the case in which the fermion has a non-zero 
mass. In that case, the resonances of the positive energy (particle) scattering solutions can 
migrate to the imaginary axis in the /c-plane and become true bound states, just as in the 
non-relativistic case. As the potential is deepened still more, however, the bound states 
eventually turn around and re-enter the continuum, this time as resonances for the negative 



energy (antiparticle) scattering solutions |]T4|. The presence of resonances in the negative 



energy scattering solutions signals that the ground state has become charged. The situation 
is similar in the massless case, except that in this case the positive energy resonances move 
to the origin and pass immediately to the negative energy continuum without ever making 
true bound states. 

In our simple model, it is possible to calculate precisely the critical values of the potential 
at which the charge of the ground state changes discontinuously. The solutions of the "zero- 
energy resonances" of the system are normalizable solutions of the Weyl equation.[^ Setting 
ui to zero in Eqs. ( ^.6| ) and (|3.7|), we find that fi and gi decouple for r > R, giving 

// cx (3.36) 

for r > R. The normalizable solutions, then, are those which kill fi in the asymptotic 
region. After a few lines of algebra, one arrives at the following simple condition for the 
critical parameters at which the charge of the vacuum changes: 

jiiaR) = 0. (3.38) 

In particular, then, the vacuum goes from being uncharged to having charge equal to 2 when 
aR = TT. 

Using the relation p.38|), it is straightforward to calculate the charge of the ground state 



for any value of aR: for each I, one simply needs to count the number of zeros of ji{p) with 
p < aR. Fig. ^ shows a plot of the charge of the ground state as a function of aR. It is 
clear from the solid curve in this plot that the charge does indeed change discontinuously 
as aR is increased. The dashed curve in this plot shows the charge which one obtains in 
the large- volume limit for a system with chemical potential p = a, Eq. ( p.23| ). As aR is 
increased, we find that our exact result tends to this value. 

Let us now calculate the renormalized energy, W^^^\ for aR > vr. Let us define 

oo 

(4+) = £ (3.39) 



1=0 



where 



= -1(2/ + 2) r duo \5i{uj) + 6i{-uj) - 2{aRf6f\uj)\ (3.40) 
27r Jo L 



^(2/ + 2) Hduj [5,(^) + 5K-^)]-^(2/ + 2)/f) (3.41) 
/TT Jo /vrit 



These are the "would-be" bound states of the system which, in a vector-hke theory with a 
massive fermion, would have migrated onto the imaginary axis in the /c-plane and become bound 
states. 
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In the second line we have exphcitly separated out the term corresponding to the integral of 
the second Born term, /p'*, since it may be calculated analytically. (See Fig. ||.) We stress 
that for a given / both terms in ( |3.41| ) are finite. 

Fig. P shows two plots of the dimensionless combination Wi^^^^R as a function of / for 
two strong potentials. In each plot the solid dots correspond to aR = 10 and the open 
circles correspond to aR = 20. These plots differ qualitatively from the analogous plots 
in the case of weak potentials. In the present case, the energy contains a large "bump" 
for / < /crit (see Fig. ^ for comparison). For / > /crit, the energy has a power- law fall-off 
which is similar to the large-/ behaviour found in the weak-potential case. From the slope 
of the tail in the log-log plot, one sees that for large / the tail goes roughly as Q The 
presence of the bump in the plot of Wi ^^^ is directly related to the presence of a condensate. 
In all cases which we have studied, the value of / at which the energy suddenly drops is 
precisely the same value of I for which there is no more contribution to the vacuum charge; 

i. e. Zcrit = /max ~ OiR. 

The final step in our calculation is to perform the sum of ly/^^"* over / in order to obtain 
the renormalized energy, and then to study the behaviour of the renormalized energy as a 
function of aR. Fig. ^ contains plots of both the renormalized charge and the renormalized 
energy as functions of the dimensionless parameter aR. This figure represents the central 
result of our paper. The plots of both the charge and the energy have been normalized to the 
values which one obtains assuming a condensate with chemical potential fj, = a, Eq. (p.23|) 
and Eq. {^1^). 

For aR < vr, the charge is equal to zero and the renormalized energy is well described 
by keeping only the first non- vanishing term in the perturbative expansion [see Eq. (|3.34| )]. 
This limit is indicated by the dashed line in Fig. [^(b). For aR > n the energy is no longer 
well- described by the leading term in the perturbative expansion, and instead approaches 
the value which one expects for a large-volume system with chemical potential fi = a. 

Fig. ^ demonstrates very convincingly that the magnitude of the shift in energy of the 
vacuum does not increase exponentially as the potential becomes large, but rather crosses 
over smoothly to the value expected for a condensate of massless neutrinos. The numerical 
calculations become quite computer-intensive for large aR and we have not gone beyond 
aR = 100. In our opinion, there is no compelling reason to go to larger values of aR. 



IV. CONCLUSIONS 

In this paper we have presented an exact non-perturbative one-loop calculation of the 
neutrino vacuum energy in the presence of an external electroweak potential. In order to 
simplify the calculation, we have chosen to work with a spherical square well potential with 
depth a and radius R. 
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Note that a very precise determination (to better than one part in 10^) of the two terms 



contributing to Wi^^^ in Eq. ( 3.41| ) is required in order to see that the tail in Fig. ^(b) is actually 



converging. 
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The formal expression for the vacuum energy is ultraviolet divergent. This divergence 
is however limited to the leading term in the expansion of the vacuum energy in powers of 
aR and can be renormalized using the standard methodology. The higher order terms are 
free of ultraviolet divergences. For large values of aR, however, these terms are infrared 
divergent. The non-perturbative resummation of these terms has been the main goal of our 
paper. We refer to the resummed expression as the renormalized vacuum energy. 

The perturbative expansion of the renormalized vacuum energy is reliable for small values 
of aR. For aR > tt, the perturbative expansion breaks down and we find that the ground 
state contains a non-zero neutrino charge. The neutrino charge may be calculated exactly 
for any aR in our model. As aR increases, the charge is well-approximated by the charge 
of a neutrino condensate with chemical potential a. The onset of the charged vacuum at 
aR = TT is accompanied by a smooth cross-over in the vacuum energy. (See Fig. |^.) For 
large aR, the energy approaches the value expected for a neutrino condensate. 

It has been argued in Refs. that long range forces due to the exchange of massless 
neutrinos could destabilize a neutron star. Given that a finite density of neutrons gives rise 
to an effective electroweak potential for the low energy neutrinos we can directly compare 
our results with these claims. According to P], a dense microscopic system of neutrons with 
radius ~ 2 ■ 10~^ cm and a ~ 20 eV contains as much energy due to the exchange of 
massless neutrinos as the total mass of a neutron star of radius 10 km; i.e., W ~ 10^^ eV. 
By way of contrast, the infrared-sensitive terms in our calculation contribute W ~ —2.3 keV 
per massless neutrino species. As emphasized previously by Abada, et a/ 0, the erroneous 
results obtained in are due to an improper counting of the processes contributing to 
the vacuum energy. When properly resummed, these processes lead to small effects, even 
for massless neutrinos, and do not endanger a neutron star. 



NOTE ADDED IN PROOF 

In this note we address some issues raised in two manuscripts [^,25] which appeared 
after the completion of the present work. 



In Ref. [Q, Fischbach and Woodahl have questioned the validity of the effective field 
theory approach used in the present paper. Their criticism is based on a discrepancy between 
our estimate of the two-body contribution to the vacuum energy of a neutron star [see the 
discussion below Eq. ( p.32|) ], 

_ a^R^A^ (4.1) 
and their estimate [Eqs. (22) and (23) in p4| ]. 



ly(^) ~ a^R^A^. (4.2) 

It is however straightforward to reconcile these estimates. The discrepancy may be 
traced to the fact that we have used a smooth potential to model the star. This mean field 
approximation, while very well-suited for studying the infrared behavior of the system, is a 
poor approximation to a realistic star in the ultraviolet, where the coarse grained structure 
becomes important. It is easy to verify that upon introducing fluctuations on scales of order 
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1/A and repeating the analysis of Section III.C, one obtains an estimate for W^"^^ which is 
in quahtative agreement with Eq. ( [4.2[ ). 

We reiterate, however, that the higher-point contributions to the vacuum energy [i.e., 
the terms with A; > 4 in Eq. ( p^.lO )] are insensitive to the ultraviolet features of the star. As 



discussed at length in the text, these contributions must be resummed, a task for which the 
effective field theory approach in the mean field approximation is correct and perfectly suited. 
The resummed result, W^'^~^\ tends to the value expected for a condensate of neutrinos with 
chemical potential a as aR gets large. This quantity gives the leading contribution to the 
vacuum energy for a smooth potential. For a realistic star, the dominant contribution to 
the vacuum energy is likely to be due to the two-point contribution. 

In Ref. Abada, Pene and Rodriguez-Quintero have independently analyzed the 

behavior of massless neutrinos in a three-dimensional electroweak potential and have, like 
us, concluded that the exchange of massless neutrinos does not threaten the stability of a 
neutron star. However, their estimate of the vacuum energy, 

W^aR^A^ (4.3) 



differs both from our two-body contribution ( [4.1|) , valid for a smooth potential, and from 
our estimate of the contribution due to the higher-point diagrams, 

^ = (4.4) 



We believe that the estimate in Eq. ([4.3|) is incomplete. Note that this expression is linear in 
the external potential a. Changing the sign of the external potential is, however, equivalent 
to considering an antineutron star (which attracts antineutrinos) instead of a neutron star 
(which attracts neutrinos). This transformation leaves the vacuum energy unchanged. That 
the vacuum energy is an even function of a is indeed manifest to all orders in the Schwinger- 
Dyson expansion. In particular, the tadpole diagram [which behaves superficially like the 
expression in Eq. ( ^.3| )] vanishes identically when both the positive and negative energy 
eigenstates are taken into account in the summation over modes. The first non-vanishing 
contribution is then the quadratic term as discussed above. 
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FIG. 1. Plots of the phase shift for aR = 1 and I = 0. (a) The upper (lower) solid curves 
give the particle (antiparticle) phase shifts, while the dashed curves give the sums of the first and 
second Born terms, (b) This curve shows the sum of the particle and antiparticle phase shifts with 
twice the second Born term subtracted. 
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FIG. 2. Plot of -{21 + 2)/^^''^ vs. I. 
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FIG. 4. Plot of the charge as a function of aR. The solid curve gives the exact charge, which has 
periodic jumps, and the dashed curve gives the charge expected for a condensate in the large- volume 
limit. 
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FIG. 5. Plot of II ' vs. I. 
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FIG. 6. Plots of -W/^+^i? vs. I for aR = 10 (solid dots) and aR = 20 (open circles). The 
slopes of the tails in the log-log plots are approximately —2 in each case. 



26 



1.0 
0.8 
I 0.6 - 
0.4 
0.2 

0.0 L 

10 



10" 10' 
aR 



(a) 



10^ 



1.0 
0.8 
10.6 
i 0.4 
0.2 
0.0 



10" 



(b) 



10" 



10' 



10^ 



aR 



FIG. 7. Plots of (a) the charge and (b) the energy as a function of aR. The dots give the results 
of our exact (non-perturbative) calculations. Both the charge and the energy are normalized to 
the contribution expected from a condensate in the large- volume limit. The lower dashed curve in 
(b) gives the value obtained in the limit of small aR (see Eq. ( 3.34| )). 
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